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We discuss the indirect and direct role of the short-lived resonances as probes of QGP freeze- 
out process. The indirect effect is the distortion of stable single particle yields and spectra by 
contributions of decaying resonances, which alter significantly the parameters obtained in fits to 
experimental data. We than discuss the direct observation of short-lived resonances as a probe of 
post-hadronization dynamics allowing to distinguish between different hadronization models. 

PACS numbers: PACS: 24.10.Pa,12.38.Mh,25.75.-q 



o 
o 

(N 

o 
O 
in 



> 
o 
o 

o 



X 



I. INTRODUCTION 

We would like to understand the dynamics of quark- 
gluon plasma (QGP) phase breakup into individual 
hadrons and show here how this can be done using 
hadron resonances. To convert QGP into hadrons we 
employ the Fermi statistical model of particle produc- 
tion 01 S IS S IE ■ This approach has been used ex- 
tensively in the field of relativistic heavy ion collisions. 
Particle abundances and spectra both at Super Proton 
Synchrotron (SPS) [1 H 11 E lil|Jl3Jlland Relativis- 
tic Heavy Ion ColHder (RHIC) glElil 113, E, 13 
energies have been analyzed in this way. In all this stud- 
ies one relies on presence of a distinct condition in tem- 
perature at which the rate of particle interactions goes 
to zero (freeze-out). The HBT _2Qj measurements and 
the success of statistical models suggest that the freeze- 
out process at SPS and RHIC energies is indeed much 
faster than previously expected, being close to the explo- 
sive hadronization limit of instantaneous emission and 
negligible post-emission interactions . 

The obvious way to test this model experimentally is 
to obtain experimental measure of the time lapse be- 
tween hadronization and end of hadron-hadron interac- 
tion. While HBT has always been considered the ideal 
probe to do it, the fact that hydrodynamic models at 
present fail to describe HBT data 12^ , together with the 
problems associated with emission from an interacting 
gas with non-zero mean free path |23j | make the search 
for a different probe of freeze-out time necessary. Since 
hadronization is a fast process, resonance decay products 
have an appreciable chance of escaping without rescatter- 
ing and thus resonance yields can be measured directly, 
and in an analysis of these results one can obtain infor- 
mation about hadronization dynamics. 

Therefore we consider direct detection of hadronic 
short-lived resonances as an alternative probe of freeze- 
out dynamics. We find that the short-lived resonances, 
detectable through invariant mass reconstruction ,2^ ^E 
|2^ |27| are natural candidates for freeze-out diagnostics 
since their lifetime is comparable to the hadronization 
timescale and the lifetime of the interacting HG. Res- 
onances usually have the same quark numbers as light 
particles, making their yield compared to the light par- 
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FIG. 1: (Color online) Relative resonances contribution to 
individual stable hadrons for three particle freeze-out tem- 
peratures. 



tide independent of chemical potential. The rich variety 
of detected resonances 25] includes particles with very 
different masses and widths, allowing us to probe both 
production temperatures and interaction lifetimes in de- 
tail. Fig. n shows what percentage of observed light par- 
ticles comes from the decays of heavier resonances (quite 
a few of them experimentally observable). As can be seen 
in figure n this resonance contribution is significant, and 
varies appreciably with both particle type and tempera- 
ture. 



II. RESONANCE INFLUENCE ON PARTICLE 
SPECTRA 



Statistical hadronization model assumes that, in the 
rest-frame with respect to collective matter flow, parti- 
cle spectra are given by the (entropy maximizing) Boltz- 
mann spectrum. In addition, the phase space occupancy 
by particles is affected by the space-time geometry of the 
particle emission surface. The particle spectrum in the 
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laboratory rest frame is given by, 




niT coshye '^^^^ "^"^ , (1) 
(2) 



Taking feed-down from resonances into account can 
be tedious numerical task. To simplify the situation we 
assume that, in a decay of the form. 



(3) 



dynamical effects in the decay average out over a statis- 
tical sample of many resonances. In other words, in the 
rest frame comoving with the "average" resonance, the 
distribution of the decay products will be isotropic. If 
more than two body decays are considered, this calcula- 
tion becomes more involved . For the general N-body 
case, evaluation is better left to Monte Carlo methods 

The rate of particles of type 1 (as in eq. produced 
with momentum pi in the frame at rest w.r.t. the reso- 
nance will then be given by the Lorenz invariant phase 
space factor of a particle of mass Mi and momentum pi 
within a system with center of mass energy equal to the 
resonance mass M^, 
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Here h is the branching ratio of the considered de- 
cay channel. To obtain laboratory spectrum we need 
to change coordinates from the resonance's rest frame 
{p* ,E*) to the lab frame (p, E). 

In the case of the 2-body decay p* , E* are fixed by the 
masses of the decay products. 
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Putting the constraints in eq. 
after some algebra 113 



(O into eq. Q one gets, 



dN 



drri^-^dyi Anpl Jy 
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Y+ i.Mt+ j2 f.T 

dvJ dM^,J '^^^ 
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(7) 

^P'^rP'ti - {MrE*j^ - MTRmri cosh Ar)2 ' 
l\Y = YR-yi. (9) 

J is the Jacobian of the transformation from the reso- 
nance rest frame to the lab frame, and the limits of the 
kincmatically allowed integration region are: 
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FIG. 2: profile for a fit to WA97 Ks, A, A, H, H, Vl, Q. data 
at 158 GeV Isi]. Solid line: resonance decays included. Dot- 
ted line: resonances reequilibrated. 
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This calculation assumes resonances are produced 
through the same statistical model which is fitted to the 
momenta of the daughter particles, and their decay prod- 
ucts reach the detector with no further interaction. How 
the data is described in this approach is seen in figure 
121 for the most peripheral, and the most central reac- 
tion bin. Solid lines show the no-rescattering fit qual- 
ity for particle spectra. When rescattering of the decay 
products is important the shape of the spectra are as if 
there were no resonance decays. This case is presented 
by dashed lines in figure |2] Clearly there is a much more 
physically significant spectra fit possible if we allow that 
hadron resonances are present and their decay products 
do not rescatter. 

Inclusion of resonances in the spectral fit requires no 
extra degrees of freedom, and the improvement of 
of the fit of the SPS hyperon spectra seen in figure [21 
|3l| convincingly favors models in which resonances are 
present within the fireball in the abundances predicted 
by a single freeze-out temperature model. Not including 
spectral contributions of resonances raises the by a 
factor of 3 and removes the minimum's significance. 



3 



We have not yet found a particle type for which consider- 
ation of resonances does not lower the x^- Moreover, we 
find considering possible mass or width shift that such 
effects significantly increase the "ill , and thus there is 
another evidence that rescattering, which usually shifts 
masses and widths 1331, is minimal. 



III. RESONANCE YIELDS AS PROBE OF 
HADRONIZATION DYNAMICS 

The fact that resonances have been found to give the 
contribution to particle spectra predicted by the statisti- 
cal model has motivated the direct search of resonances 
through invariant mass reconstruction p^ . I25I |27| . It 
has been found that while statistical models are able to fit 
the K* / K ratio within error j^, the A(1520) is strongly 
suppressed at both SPS and RHIC energies [2^. This 
puzzling result suggests that both emission temperature 
and re-scattering might make a significant contribution 
to resonance abundance, and their roles need to be dis- 
entangles. This is possible by considering abundances of 
two resonances with different masses (which constrains 
the temperature) as well as lifetimes (which probes the 
role of non-equilibrium in-medium effects, i.e. rescatter- 
ing). 

To obtain a quantitative estimate, we have calculated 
the ratios of {K* +1{*)/Ks and A(1520)/A using the 
statistical model. In both cases chemical potential cor- 
rections are negligible since the particle's chemical com- 
position is the same and thus we have: 



N* 



n{m*, T) 



N + N* n{m*,T) + n{m,T)'' 



n{m, T) = m TK2 { — 



(10) 



(11) 



We then evolved in time the ratios using a model which 
combines an average rescattering cross-section with di- 
lution due to a constant collective expansion. In this 
model, the initial resonances decay with width T through 
the process A^* D. Their decay products D {D{t = 
0) — 0) then undergo rescattering at a rate proportional 
to the medium's density as well as the average rescatter- 
ing rate. The final evolution equations then are 



dN* 
IT 
dD 



= -TN* 



R, 
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.(13) 



V is the expansion velocity, Ro is the hadronization ra- 
dius, pj — Tij (mj , T) the initial hadron gas particle den- 
sity and {(JDjVjjj) is the particle specific average flow and 
interaction cross-section. 

In this calculation, we have neglected the regeneration 
term R oc {(jj^^^VDi)pi, since detectable regenerated 
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FIG. 3: Observable relative resonance yields as a function of 
temperature for a given length in time of interacting hadron 
gas phase. Horizontal and vertical lines give experimental 
results. 



resonances need to be real (close to mass-shell) particles. 
Figure 13 shows how the ratios of S*/A and K* /K evolve 
with varying hadronization time within this model. 

It therefore becomes apparent that measuring two such 
ratios simultaneously gives both the hadronization tem- 
perature and the time during which rescattering is a sig- 
nificant effect. Figure ^ shows the application of this 
method to if*, A(1520) and S*. 

The special role of A(1520) suppression is evident. 
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FIG. 5: Spectral ratios as function of m± for an array of 
freeze-out surface conditions. 



measured E* abundance. 



IV. RESONANCE SPECTRA AS DIRECT 
PROBE OF FREEZE-OUT DYNAMICS 



FIG. 4: Projected combined relative yield diagrams over a 
mesh of temperature and lifespan. Horizontal and vertical 
lines give experimental results. 



A(1520) is a very peculiar resonance, since unlike the K* 
and S* it's extra spin is believed to originate from inter- 
quark orbital angular momentum (L=2). It is therefore 
particularly susceptible to in-medium effects which sup- 
press it's yield or enhance it's width |33,I3- this is 
the case, our model is able to account for existing ob- 
servational data and makes definite predictions for the 



The near-independence of resonances on chemical po- 
tential makes their direct detection a powerful tool for 
examining further aspects of freeze-out dynamics. It is 
apparent from the previous section that the ratio of the 
resonance to the daughter particle with the same number 
of quarks is extremely sensitive to freeze-out tempera- 
ture. If the m_L dependence of this ratio can be observed 
further freeze-out parameters can be extracted [s^. For 
a purely thermal source, in which 
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the observed resonance ratio should be a step function, 
with N* /N — for m± < m* and the ratio of degenera- 
cies for the resonance g* / g for m± > m* . 

However, a non-trivial flow profile and emission geom- 
etry can signiflcantly change this dependency. For this 
reason, this ratio is an extremely sensitive probe of both 
flow profile (the transverse flow as function of the ra- 
dius) and hadronization dynamics (how emission time 
varies with space, parametrized by w/ in eq. Figure 
|S1 shows this ratio (calculated by dividing two expres- 
sions of the form of eq. but different masses) done 
for 77' — > 77/ry — > 77 and {K* + K*)/Ks- It is appar- 
ent that this measurement is indeed a sensitive probe 
of freeze-out dynamics and flow. Moreover, flow effects 
are separated from freeze-out dynamics, something which 
"normal" hadronic spectra can not do effectively 36]. 

In conclusion, we have presented here an outline of 
the use of resonances as diagnostic tools in the study of 



freeze-out dynamics. We have shown that statistical flts 
to particle spectra require an admixture of resonances 
consistent with thermal predictions, which strongly sug- 
gests negligible post-hadronization rescattering dynam- 
ics. We have described how the study of resonances yields 
and spectra can yield the hadronization temperature, the 
timescale of thermal freeze-out, and the freeze-out geom- 
etry dynamics. We expect these models to be constrained 
and developed further as resonance experimental data be- 
comes more precise. 
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